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In this paper one examine analytical solutions for flat and non-flat universes composed 
by four components namely hot matter (ultra-relativistic), warm matter (relativistic), cold 
matter (non-relativistic) and cosmological constant. The warm matter is treated as a reduced 
relativistic gas and the other three components are treated in the usual way. The solutions 
achieved contains one, two or three components of which one component is of warm matter 
type. A solution involving all the four components was not found. 
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I. INTRODUCTION 

At recent years cosmology has lived a gold age. Many observational techniques are being 
developed and they are producing a lot of data about our universe. We can cite for example 
techniques involving observation of supernovas [IHS], detection of cosmic microwave background 
anisotropies [6] and surveys of galaxies [7H9]. These new data have generated a picture where the 
universe is constituted by five components: radiation, neutrinos, baryonic matter, dark matter and 
dark energy. In the simplest models the components evolve independently which means each one 
has its one equation of state. 

In the standard approach the evolution of the universe is divided in eras, each one dominated by 
one component. Separate solutions are used for each era and intermediate periods are connected 
matching the initial conditions. This approach, although approximated, describes the history of the 
universe in a simple way. However, the precision of new data, make it desirable to have solutions 
as complete as possible. 

Analytic solutions involving cold matter, radiation and cosmological constant were studied in 
|lUH13j . Nevertheless, none of these articles consider relativistic components in their analysis. In 
principle, if we look for solutions that take into account relativistic particles it would be necessary 
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to deal with equations of state containing modified Bessel functions. This would make it very 
difficult to obtain analytical solutions. Fortunately, in 2005 it was proposed a simpler formulation 
to describe relativistic particles [Hj. This formulation, known as reduced relativistic gas (RRG), 
is able to represent a gas of relativistic particles with good accuracy. Besides, the RRG model is 
simple enough which allows us to search for analytic solutions to Friedmann equations. 

The first analytic solutions containing the RRG were found in |14j . Here one continue this 
work extending the analysis for cases involving RRG with other components. It is discussed 
solutions for a universe composed by a RRG component plus non-relativistic matter, radiation 
and/or cosmological constant. These kind of solutions are important whenever we want to describe 
a universe which has a component with relativistic behavior. Good examples are models which 
involves warm dark matter |15} H6] . 

The paper is organized as follows. In section 2 is given a summary about RRG focusing in its 
connection with standard cosmology. The content of this section is a resume of sections 2 and 3 
of [2]. In section 3 is presented analytic solutions for a universe containing one, two and three 
components where one of these components is modeled by RRG. In general, the flat and non-flat 
cases were treated separately. The final comments and further perspectives are given in section 4. 



II. REDUCED RELATIVISTIC GAS IN THE STANDARD COSMOLOGY 



The RRG is a simple model for a relativistic ideal gas of massive particles and it was first 
introduced in |14j . The idea behind of this model is to use the kinetic theory of gases attached 
with relativistic concepts. Using standard considerations which relate the transferred relativistic 
moment by particles with the pressure p produced in a wall, allow us to write 



mv 



P 



3V -v 2 /c 2 

where V is the volume, m is the mass and v is the relativistic velocity. 

Supposing that all particles have the same relativistic kinetic energy e one can rewrite the 
equation above as 

„2\ 21 



P 



1 



nmc 



(1) 



where p = ne is the energy density and n is the numerical density of particles. Note that pd is the 
energy density of non-relativistic particles and thus it is proportional to V, i.e. pd = p\V . It is 
easy to see that if p ~ pd (non-relativistic particles) the equation of state ([TJ reduces to p ~ 0, and 
if p 3> pd (ultrarelativistic particles) the equation of state (EoS) becomes p = p/3. 
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It is instructive compare ((!]) with the correct EoS derived by the statistical mechanics of ensem- 
bles. Computing the partition function of classical relativistic ideal gas we can determine p and p 
as functions of n and kT: 

p = nkT 

2 K 3 {mc 2 /kT) 
K2 (mc z /kl ) 

where K v is a modified Bessel function of index v. Combining this two equations we obtain 

K2 {Pd/P) 

At first sight, ([I]) and ([2]) are completely different. However, a numerical comparison between 
these two equations was performed in p3] and there it was shown they are quite similar. Indeed, 
the difference between and ([2]) is at most 2.5% and becomes negligible at ultrarelativistic and 
non-relativistic regimes. Thus, is a good approximation for the EoS of classical relativistic ideal 
gas with the great benefit of being much simpler than ^ . 

In order to use ([I]) as an EoS for a relativistic component of cosmic fluid it is necessary deter- 
mined how the energy density depends on scale factor. This is performed writing the covariant 
conservation law in terms of volume 

< | = - (/ ^ ) where V ~ a 3 , (3) 
and replacing ([!]) in Q. Solving the differential equation we obtain 

where the initial condition used was prrg («o) = \/ Pi + p\ ■ 

Analyzing the last equation we can associate the constants p\ and P2 as the energy densities 
of dust and radiation respectively. Indeed, if we take p2 = 0, Q scaling as a dust-like component 
which means p ~ a~ 3 . And if we take p\ = 0, Q scaling as a radiation-like component which 
means p ~ a -4 . Although, Q reproduces these two behaviors it is qualitatively and quantitatively 
different from a cosmic fluid composed by dust and radiation. In the first case, we have a single 
relativistic component represented by Q, and in the second case, we have two distinct components 
whose the energy density is given by 



p m+1 (a) = pt (— ) +P2["^) ■ (5) 
A numerical confrontation between Q and ^ is shown in figure [TJ 
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FIG. 1: The energy density of prrg (full line) and p m + 7 (dashed line) in terms of a/ag- In both case it was 
chosen p\ = p2 with the initial condition p(ao) = 1. Observe that the transition among the radiation and 
dust behavior is smoother for prrg- 

III. ANALYTIC SOLUTIONS 



Suppose that the cosmic fluid are composed by four independent components namely radiation 
(7), cold matter (CM), cosmological constant (A) and warm matter (WM). Thus, the first 
Friedmann equation in units of c = 1 results in 

2 



8vrG 



P70 



+ PC MO 



«0 



+ 



A 



+ 



K 

^2' 



(6) 



a J \ a J 8irG 

where k is the spatial curvature and the WM is modeled by the RRG. Usually, the component 
of WM represents a warm dark matter (as proposed in |15| [16]). Nevertheless, another physical 
possibility is use the WM component to describe neutrinos. 

The main goal of this work is to study the analytical solutions linked with Q. To perform this 
it is convenient to define the following quantities: 

Pcmo ~ A 



fi 7 o 
fii 



P70 

Pc ' 

Pi 

Pc 



pi 



Pc 

and O k0 



—k with p c 
fin 



3H$ 
8vrG' 



(7) 
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Since ^ is a separable differential equation of first order it can be written as an integral in the 
scale factor. Thus, using ([7]) we obtain, 

If ada , . 

t(a) = ±— (8) 

U ° J - 'n i0 + QcMoa + Q A a 4 + O x [a 2 + 6 2 ] 1/2 - Q K0 a 2 



where ao = 1. The sign will be chosen so to get always an expanding universe. 

The approach adopted here is to search solutions involving one, two and three components 
with and without curvature 1 . The integral in ^ can be written as a sum of integrals whose its 
structures are of type 

"f (a, ^/P(a)) 

where P (a) is a polynomial. Associated with this integral we have three possibilities [T7] : 

(i) If P (a) is at most second-degree polynomial and F has a simple structure then the solution 
of Q can be expressed in terms of algebraic functions. In this case, an explicit expression 
sometimes can be found. 

(ii) If P (a) is third- or fourth-degree polynomial and/or F has not a simple structure then the 
solution of ^ can be written in terms of elliptical integrals (see appendix [AJ . In this case, 
the solution is only implicit. 

(iii) If P (a) is more than fourth-degree polynomial then it is not possible to obtain a solution 
for ^ - e.g. the integral with all four components. 

Before to move on for the specific cases it is noteworthy that the solutions without WM is not 
treated in this paper. This kind of solution was extensively studied in [10H13] . 

A. Solutions with one components 

The first and most simple case is when only one component (WM) is present. Thus, ^ is 
reduced to 

If ada 

t(a) 



1 An analytic solution involving all the four components was not achieved. 
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The solution for flat curvature (k = 0) is given by 



t f (a) 



2 r -> ,21 3/4 



[a 2 + b 2 Y^ + i^a f (t) 



ttiH (t-t) 



4/3 



b 2 



(9) 



where i is the integration constant and fii (l + b 2 ) 1 ^ 2 = 1. This result was first derived in [14J and 
it is presented here only for completeness. 

The structure of solution for non-flat cases (k = ±1) is as complicated as the solutions involving 



7 and WM. Therefore, it will be presented in the section IIIB 



B. Solutions with two components 



Solutions with two components are of type (WM, 7), (WM, CM) and (WM, A) with and with- 
out curvature. Because of their complicated structure, the cases (WM, A) with and without cur- 



vature will be treated in the section IIIC Let's perform the analysis of the two other cases. 



1. Universe with WM and 7 

For an universe constituted by WM and 7, the equation ^ is reduced to 

ada 



t{a) = ± wJ 



ft 7 o + fti [a 2 + b 2 ] 1/2 -n K0 a 2 



The solution for flat curvature is given by 



tf a) 



Ha 



2Va 2 + b 2 - 4fi 70 ) Jn^ + QWa 2 + b 2 



3^i 



+ t 



(10) 



where t is the integration constant and fi 7 o + f2i (l + & 2 ) = 1- Note that if we take S7 7 o = the 
result Q is recovered. This results was first derived in jT3j and again it is presented here only 



for completeness. Unfortunately, equation (10) can not be inverted and thus it is not possible to 
derive an explicit solution. 

For positive curvature the solution is 

-1 



t p (a) 



|^ko| 3//2 Hq 



fii -2|fi«o| Va 2 + b 2 
2y P 



+ t 



(11) 



where 



Vv 



^l Q a 2 + Oi |fi K0 | V a 2 + b 2 + fi 70 |O w0 | 
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And for negative curvature the solution is 

1 



t n (a) 



I ^kO 1 3 ^ 2 Hq 



y n - — In 



fii + 2 |n« | \/a 2 + & 2 + 2y r , 



+ t 



(12) 



where 



y n = \/ tt 2 K0 a 2 + Oi |S1 k0 | \/a 2 + b 2 + 7 o |n« |- 

_ 1/2 

In both cases t is an integration constant and Sl 7 o + ^i (l + b 2 ) — S7 k q = 1- As it should be, the 



equations (11) and (12) are reduced to (10) in the limit \Cl K o\ — > 0. This statement can be verified 



expanding the functions Arctg [...] and In [...] in powers of ^/|f2 re o| until third order. 



2. Universe with WM and CM 



For a flat universe composed by WM and CM, the equation ([8]) becomes 

1 f ada 



Using the definitions 



y 2 ^ 



t f (a) 



rb 



QcMoa + fi! [a 2 + b 2 } 



1/2 



a + Va 2 + b 2 ' 



= Arc sin (y) and r = 



CMO 



Ox + O 



CA/0 



the solution for ( 13 ) is written as 



t f (a) 



V2P 



(Oi + VLcmo) 



1/4 



;i-r 4 )F («/>,-!) + 



(r 4 



(13) 



+ t 



(14) 



where t is the integration constant and F ((f), m) is the elliptical integral of the first kind (see 

1 /2 

appendix A). Besides, we have the following constraint CIcmq + ^i (l + b 2 ) = 1. This solution is 
valid for all physical values of Qcmo and Oi avoiding Qcmo = ^i- If we take ^cmo = then r = i. 
Using this result and performing some simple algebra we recover the solution ^ for a universe 
composed only by RRG. 

The structure of solution for non-flat cases is as complicated as the solutions involving 7, WM 
and CM, and thus, they will be discussed in the next section. 



C. Solutions with three components 



Solutions with three components are of type (WM, CM, 7), (WM, 7, A) and (WM, CM, A). It 
is not possible to achieved an analytic solution for the case (WM, CM, A) because the term inside 
the square root in ([8]) is a polynomial of degree greater than four. Let's perform the analysis of the 
two other cases. 



s 



1. Universe with WM, CM and 7 

Suppose an universe composed by WM, CM and 7. In this case, equation ^ becomes 

t(a) = ±— 1 (15) 
U ° J ^n j0 + n C M a + fix [a 2 + 6 2 ] 1/2 - O K0 a 2 

This integration can be solved for flat and non-flat cases, but for k = ±1 the expressions are 
rather complicated. For negative and positive curvatures the solutions involve four and seven 
elliptic integrals respectively. Besides, each solution has constraint related with the cosmological 
parameters. Thus, they will not be presented in this paper. 

On the other hand, the flat case is relatively simpler than non-flat cases. Indeed, if we make 
the following definitions 



b (Oi + n C Mo) b(fli + ttcMo) 

2 r b a ■ / \ i s 

y = , <ft = Arc sm (y) and m = -, 

a + \Ja 2 + b 2 r 



the k = solution of (15) is given by 



1 / 2&3 r 3 j ^( y 2 _ 1) (my* _ J) 

*/ (a) ~ 



"6if V (^1 + fi CMo) I y 



■u 4 

l + 2(l+m)y 2 + ^4 



+ 



+ ( ^p- ) [(m + 2) F (0, m) - 2 (m + 1) E (0, m)] }> + f (16) 



where i is the integration constant and F ((f), m) and E (eft, m) are the elliptical integral of the first 
and second kind respectively (see appendix |A[). This solution is valid for almost 2 all physical values 
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of f2 7 o, CIcmo, ^1 and b satisfying the constraint fi 7 o + Dcmo + f^i (l + 6 2 ) =1. Nevertheless, 
the choice of sign depends on relation between FLcmo arid fii. If fipMO > ^1 (^cmo < ^1) the 
sign plus (minus) must be used. 



As it should be, the solution (16) contains the previous case involving only WM and CM. 
Taking J7 7 o = we get m = — 1 and after some straightforward algebra we recover the solution 
pl). 



2 The values b — or Qcmo = fii are not allowed. 



2. Universe with WM , 7 and A 



For an universe constituted by WM, 7 and A the equation ^ is given by 

ada 



Ha) = ± iJ 



'OyO + ft A a 4 + fii [a 2 + b 2 ] 1/2 - n K0 a 2 
It is convenient to change the variable of integration a using the relation a 2 = x 2 — b 2 . Thus, 

1 f xdx 



t(a) = ±- 



Ho^ThiJ Vx 4 - Lx 2 + Mx + P 



where 



L = 2b + 



n 



M 



fin 



and P 



n. 



i, f"'"'/ 



' ^A ^A ^A 

The new integral is not too simple but it can be solved through the following steps 3 : 

1. Rewrite the fourth-degree polynomial in terms of the roots rf 

x 4 — Lx 2 + Mx + P = (x — 7*1) (x — r 2 ) (x — r%) (x — r^) 

2. Introduce some convenient new constants: 

S = \/2(L 2 + 12P) 
W 

U = -4(L 2 + 12P) 3 + (-2L 3 + 27M 2 + 72LP) 2 



-2L 3 + TIM 2 + 12LP 



V 



s 



+ 



3(W + VE7)V3 
3. Express the roots in terms of these constants: 



(w + Vu) 1 / 3 

3^2 



+ 



2L 



; r 2 



V 



2L-V + 



2M 



W+J 2L-V 



2M 



1 








; r4= 2 





2L-V 



2M 



4. Sets the parameters m, n and the amplitude (ft as: 



r 2 - r 4 n — r 3 
n = , m = n . 



n — r4 



arcsm 



r-2 - r 3 



/ (x- 


- r 2 ) 




- n) 



These steps were first developed in [10] 
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5. The solution will then be: 



t(a) 



2(x-a) 



n 



\/(. X ~ r 2) 



a/(?"3 - r 2 ) 2 (rr - n) 



F [(A, ml + 



+ 



r 4 - r 2 



(ri - r 2 f 



\/(r± - r 2 ) 2 V (r 3 - r 2 ) - r x ) 



r II [n,0,m] 



(19) 



where x = \/a 2 + 6 2 , i is the integration constant and F (4>, m) and II (n, 0, m) are the 
elliptical integral of the first and third kind respectively (see appendix [AJ . Besides, we have 

1 /2 

the following constraint f2 7 o + Q\ + Q\ (l + 6 2 ) — f2 K o = 1- 
At this point, some features about this solution must be clarified. At first sight it seems that 



( 19 ) could be simplified. However, as (r. 



and (x — ri) could be complex numbers, any extra 



desirable simplification must be done with caution and only when the values of O 7 0) ^i and 



b are specified. Other important point is that (19) is not valid for all physical values of Sl 7 o, ^a, 
ill an d b. It happens because there is an arbitrariness in choice of which root will be n, r 2 , r% or 



r 4 . Nevertheless, the choice it was made in (18) include wide ranges for the parameters comprising 



inclusively the ACDM case. For practical purposes, a set of conditions that ensure a physical 
solution are 



0<J1a<2, 0<^i<2, 0<6<2 
7 o Ha and ^ 7 o *C 



(20) 



IV. FINAL COMMENTS 



In this paper we derived analytical solutions for a universe composed by one, two and three 
components where one of them represents warm matter. The first solution obtained is one that 
involving only warm matter. It is very simple but it serves such a guide for the complex ones. 
The next step it was to derive solutions containing warm matter plus radiation or cold matter. As 
expected, these kind of solutions are more complicated than the previous one and only implicit 
solutions were found. The most complicated solution which were achieved are ones involving warm 
matter, radiation and cosmological constant or warm matter, radiation and cold matter. These 
type of solutions, with three components, always involving elliptic integrals. Unfortunately, an 
analytic solution containing all the four components was not obtained. 
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The warm matter could mimic dark matter, neutrinos and even baryonic matter. Thus, these 
solutions can be apply in different context. For example, we can use them to analyze the effects 
of warm dark matter in structure formation [JS]. Other possibility is use them to study massive 
neutrinos in cosmology. It is noteworthy that although the results obtained concerns only to the 
background, they are also important in perturbative cosmology. Indeed, the perturbative analysis 
becomes simpler when the analytical solution for the background is known. 

Finally, it is important to emphasize that the warm matter is represented by the RRG model 
which is an approximation for a classical relativistic gas. In the context of thermodynamics, this 
approximation differs from the real situation at most 2.5% |14j . Nevertheless, none comparison 
was done at cosmological context. We expect to explore this issue in the near future. 
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Appendix A: Appendix 



Definition of elliptic integrals [T7] : 



• First kind: 



sin <j> 




• Second kind: 




• Third kind: 



sin</> 
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